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A cross-section area in in 

output voltage of a strain gauge bridge 

4 

I section moment of inertia in in 

4 

J section polar moment in inertia in in 

K gauge factor of a strain gauge, i e , unit 
resistance change per unit strain 

M bending moment in in- lb 

Q torque in in- lb 

R resistance in ohms 

T thrust load in lb 

V bridge excitation voltage 
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1 Introduction 

This project, initiated in September 1971, has as its objective 
the determination o£ propeller design parameters which yield minimum 
pro’peller noise, given constraints on propeller performance 

The study addresses both the theoretical and the experimental 
aspects of the problem Theoretical efforts have been directed at 
applying variational techniques to the noise equations Experimental 
efforts have been concerned with measurement of the propeller performance ' 
and noise of a Lockheed Y0-3A aircraft 

The following sections presents a report of the current state of 
this study 
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2 Theoretical Studies 

The objective of the theoretical aspect of the program is to deter- 
mine the optimum blade loading to minimize noise level, given the thrust, 
diameter, rotational speed and forward velocity. Both an aerodynamic model 
and an acoustic model for the propeller are required for this analysis 

2. 1 Aerodynamic Model 

The aerodynamic model at present consists of a vortex system util- 
izing a lifting line and a helical vortex sheet for each blade The blade- 
wise circulation distribution F (y) is the sought for quantity in the noise 
minimization problem 

A computer program has been written for the aerodynamic model based, 
essentially, on the analysis of Moriya (Ref. 21). The program, as it now stands 
assumes constant blade^vise section aerodynamic characteristics. Modification 
to allow variable section aerodynamics is a relatively easy task and is now 
in progress The results of the aerodynamic analysis express the propeller 
performance in terms of the Fourier coefficients for the circulation distri- 
bution r(y) A draft of this analysis is enclosed as Appendix A. Some 
further shakedown and check out of the program is required before it becomes 
fully operational 

2. 2 Acoustic Model 

The acoustic model is a simplification of that represented by the 
Ffowcs-Williams-Hawkings equation (Ref 2.2, Eq 3-16) Presently only the 
force terms are retained, resulting in a model that is essentially Gutin's 
model with radially distributed line dipoles representing the blades 
Following successful treatment of this model, volume displacement effects 


will be included. 
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The sound intensity may be expressed in normalized form, as 


I* 


J 

o o 


F ( r* , rp Y Cr*) Y ( r| ) dr-^dr* 


where y is related to the circulation and F is a known function The con- 
straint of constant thrust is expressed as 
1 

y (r*) (l-r*)dr* = 1 

0 


Applying- variational analysis to these equations provides a necessary condi- 
tion for to be a minimum 

h 

^ F(r*,rpyCrpdrJ = - ^ r^(l-r*) 

0 

where b is a Lagrange multiplier. 

Initial attempts at numerical inversion "of this integral equation 
utilized both a Legendre Polynomial expansion and a Chebyschev polynomial 
expansion for yC^*), coupled with a Gaussian integration procedure, reducing 
the integral equation to a matrix equation 

AL = C 

where L represents the unknown values of y at selected points on the blade 
Numerical inversion of this matrix equation was not successful In all cases 
tried, the coefficient matrix A was ill-conditioned 

The current approach involves a Fourier transform of the integral 
equation, implying the assumption that the noise is periodic in the blade 
frequency, with a Doppler shift This is consistent with the acoustic fax- 
field assumption 

Some success is indicated with this approach, however convergence 
IS exceedingly slow and, as a consequence, required computer time is large 
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Appendix B contains some details of the earlier analysis for a 
simplified case in \tfhich only the thrust component of blade lift was included. 

Appendix G contains some details of the current analysis utilizing 
Fourier Transforms. 

2 3 References 

2 1 Moriya, T., "Selected Scientific and Technical Papers," University 
of Tokyo, Tokyo, 1959 

2.2 Goldstein, M E , "Aeroacoustics," National Aeronautics and Space 
Administration, Washington, D C , 1974 
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3. Experimental Activities 

All experimental efforts have been directed toward in-f light measure- 
ments of propeller torque and thrust and fly-by sound levels, utilizing a 
Y0-3A aircraft 

3 1 Aircraft Preparation 

Since the Y0-3A aircraft was received in disassembled, worn and 
damaged condition, a good deal of time and effort has been spent rendering 
It airworthy. Much of this effort was expended in obtaining maintenance 
and operating manuals, shop drawings, etc. The cooperation in these endeavors 
of Mr Harold Schuetz, USAAVSCOM, St Louis, and Mr. David Schnebly, Lockheed, 
Sunnyvale, is very gratefully acknowledged 

Major repairs have included a top overhaul on the engine, repair of 
a fractured vertical fin-fuselage attachment fitting, repair of aileron 
and wing tip damage, replacement of elevator attachment pins, manufacture 
of special Jigs for checking drive shaft/prop shaft alignment and belt ten- 
sion, modification of the forward instrument panel, installation of VHP radio, 
removal and repair of the muffler system, and other items Photographs of 
the aircraft and the propeller drive systems are shown in Figures 3 1 and 
3 2 

3.2 Measurement of Propeller Performance 
3 2.1 Introduction 

An essential element in the development of a low-noise propeller is 
the accurate measurement of its propulsive performance This means that the 
thrust available, the thrust power available, and the shaft power required 
must be measured, preferably in flight For the subject studies, a 
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Lockheed YO-3A aircraft is available for use as a test platform A signifi- 
cant part of the program has been devoted to instrumenting this aircraft 
for the in-flight measurement of thrust and power 

Two approaches to these measurements were considered in this study 

a) Airflow Measurement 

By pressure surveys of the flow field behind the propeller, 
determinations of the thrust and torque acting on the propeller 
are possible (Ref 3.1) However, this method requires correc- 
tion for the interference effects of the fuselage The accuracy 
of this method was judged insufficient for the purposes of this 
proj ect 

b) Direct Thrust and Torque Measurements 

The YO-3A’s propeller is mounted on a cantilevered pro- 
peller shaft (Figures 3 2 and 3 3) This shaft is driven 
through a speed reduction drive system employing pulleys and 
belts Because of this drive system, the engine is isolated 
from the shaft and load cell measurements on the engine mounts 
would not yield propeller thrust data Such measurements could 
produce approximate torque measurements if pulley, belt and 
propeller shaft bearing losses could be determined 

As shown in Figure 5. 3, approximately 9 -inches of the pro- 
peller shaft IS easily accessible The strains induced in this 
shaft are direct measures of the thrust and torque transmitted 
to the propeller Thus, direct instrumentation of the propeller 
shaft (say, with strain gauges) can, in principle, yield the 
desired information This approach was selected 
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There is one major difficulty involved in this approach, i e , 
the thrust related strains are much smaller than the torque re- 
lated strains and are very small in an absolute sense No mechani- 
cal modifications to the shaft (e g weakening the shaft in ten- 
sion) were attempted because of airworthiness questions involved 
Instead, it was decided to attempt to measure the shaft strains 
directly 

322 Propeller Shaft Loads 

The Y0-3A aircraft is powered by a Continental IO-360A engine rated 
at 210 hp at 2800 rpm It drives a 3-bladed, constant-speed propeller 
through a sheave and belt system that provides a 3 33 1 speed reduction 
ratio. The blade pitch is controlled by a Woodward governor that supplies 
pressurized engine oil to the propeller through the hollow propeller shaft 
Based on data given in Reference 5 2 and assuming no drive system 
losses. It IS estimated that a maximum of 186 hp is available at the pro- 
peller shaft to drive the propeller This, together with a propeller shaft 
speed of 841 rpm, yields a maximum torque to be transmitted by the shaft 
of 13,944 in -lbs. 

Again based on data from Reference 3 2, it appears that the thrust 
required in level flight, at sea-level, will be in the range from 250 to 
400 lbs. Although performance data are not available for the propeller, 
the assumption of an efficiency of 50% at the sea-level stalling speed of 
61 knots yields an estimated maximum thrust available of about 500 lbs. 

The thrust load instrumentation has been designed on the basis of a maximum 
thrust load of 500 lbs, but can measure thrust loads up to 1000 lbs 

Besides the torque and thrust loads , the shaft is subj ect to loads 
from three other sources 
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aO Bending Due to Propeller Weight 

Bending stresses due to the weight o£ the propeller are estimated 
to be due to a 1500 in -lb moment acting about the location of the 
instrumentation on the shaft 

b) Internal Oil Pressure 

The internal oil pressure produces an axial tension stress due 
to the action against the hub piston. It also produces a hoop 
stress (tension) that, through Poisson’s ratio, induces a com- 
pressive axial stress The net effect is a small axial tension stress. 
This net stress is about 28 psi for an oil pressure of 100 psi; this 
IS equivalent to an additional thrust load of about 53 lbs 

c) Centrifugal Force 

At a shaft speed of 841 rpm, the outer surfaces of the 
shaft has an acceleration of about 24 g's The resulting axial 
stress IS negligible compared to the other shaft stresses 
Table I summarizes the stresses and strains in the propeller shaft due to 
thrust, torque, and bending Both the centrifugal stresses and the net 
oil pressure stresses as calculated were negligible, however, their effects 
Will be evaluated during calibration of the instrumentation The ratios 
of the maximum strains taken from Table I are 

^T ’ ^0 . 1 6 6 . 42 4 

max max ^ax 

These ratios show the crux of the instrumentation problem, i e., the very 
low thrust strains in the presence of much larger bending and torque strains 
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323 Strain Transducer Arrangements 

The sensor choice was dictated by the need to distinguish between the 
combined loads present, as well as the need to measure them reliably and 
accurately The sensor must also function over the range of environmental 
conditions met in aircraft flight. Strain gauges were chosen as the sensors 
Most strain gauge work is done in the 50-500 microstrain range, 
within this range it is possible to measure changes of 1 or 2 microstrains 
However, the thrust strains on the shaft are 8. 85 microstrains full scale, 
see Table I. This low strain level poses problems in the thrust measure- 
ment The torque and bending strain levels are high enough to present no 
difficulties in reliable measurement by foil strain gauges 

Figure 3,4 outlines the decision process followed in choosing the 
method of measurement and selection of sensors Two main avenues of approach 
were made for the thrust measurement 

1) Intensification of the low level thrust strains by mechanical 
means, using foil strain gauges as sensors. 

2) Direct shaft strain measurement of thrust strains Foil 
strain gauge and semi-conductor strain gauge evaluations were 
made and compared The semi-conductor gauges are 75 times as sen- 
sitive as foil gauges but have a high inherent temperature effect 
on sensitivity as compared to foil gauges 

Wheatstone bridges, with an active strain gauge in each arm were 
chosen If all four arms experience equal temperature changes, and if all 
four gauges are perfectly matched then temperature effects are automatically 
cancelled. Figures 3 5 and 3.6 illustrate the directions of the strains 
produced by the thrust, bending, and torque loads and 'the bridge arrangements 

used for measuring them. Each of the strains can be measured independently 
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of the other two and, in each case, the bridge output exceeds the output of 
a single active gauge 

When the shaft is subjected to thrust, bending and torque, the in- 
duced strains act along certain directions, as shown in Figure 3.5 Thrust 
and bending loads produce tensile and compressive strains in the axial 
direction Due to the Poisson effect, axial tensile and compressive 
strains generate compressive and tensile lateral strains. The shear strains 
induced by torque cause tensile/ compressive strains on mutually perpendicular 
axes inclined at 45° to the shaft axis, the sum of the components of torque 
strains in the axial and laterial directions are zero The sums of the 
thrust and bending strain components, in the torque strain directions, are 
not zero but their effects on the torque bridges will cancel 

Figure 3 6 shows the three arrangements of strain gauges used to 
measure thrust, bending and torque on the shaft, taking advantage of the 
directional properties just discussed Referring to this figure, we have 
Thrust, [2(l+u)e.j,J V 

Torque [4 Eq] V 

Bending E^ = [4eg] V 

I 

Ignoring second order effects, non-linearities and cross-sens\tivi- 
ties of the gauges, and assuming the gauges are perfectly matched and aligned 
in the required directions, each of these bridge arrangements reacts only 
to one type of loading and is insensitive to the other two loads 

3.2 4 Mechanical Intensification of Thrust Strain 

Larger thrust strains can be obtained by adding a load path parallel 
to the propellant shaft and incorporating in this added path a short, weak 


link 



Such an arrangement is shown in Figure 3. 7, The propeller shaft (1) carries 
most of the thrust load The parallel load path (made up of the support link 
C2) and the weak link (3)) carries only a small part of the thrust load, the 
amount is determined by the relative strengths of the weak link and the pro- 
peller shaft However, both load paths stretch equally under the load 
Because most of the stretching of the parallel path occurs in the short ^veak 
link, its strain is much greater than that of the shaft, i e , >> 

An analysis of this strain intensifier arrangement yields the following 


intensification ratio 




'■A F 


where 


e = strain, micro-inches per inch 
Z = length 

A = cross-sectional area 
E = modulus of elasticity 

The subscripts 1, 2, and 3 refer to the components as shown in Figure 3.7 

and is the shaft strain in the absence of the intensifier If the tensile 
0 

strength of the weak link is made much less than those of the support link, and 
the propeller shaft such that 


and 


'^1 

1^1 


^2 ^3^3 

C/) « 1 

^3 ^ 2^2 


then the intensification ratio is approximately 
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The system used to study the performance of the strain intensifier is 
shown in Figure 3.8. The system consisted of two weak link assemblies, 180° 
apart, mounted on a full-scale model of the exposed portion of the propeller 
shaft. The links were fabricated of aluminum and were each instrumented 
with 4-arm strain gauge bridges. Using the smallest gauges available, a 
link length of 0.2 inches was achieved, yielding a length ratio, of 

40. The thickness and width of the links were 0.010 and 0.400 inches, 
respectively. The tensile strength ratios were 


= 0.003 

^ 1^1 


^ 3^3 

V2 


= 0.001 


The test assembly was installed in a lathe (Figure 3.9) and could be 
loaded simultaneously in torque, bending, and thrust. Loadings up to 40% 
of full-scale were used. Typical results of the thrust calibration are 
sho\m in Figure 3.10, and the overall thrust and the torque interactions 
with thrust results are summarized in Table II. It is evident that the 
torque interaction is very large. Over the full-scale range of loads to be 
met in flight, the weak link system would have an output, due to torque 
interaction, of 860% and 212% of the full-scale thrust output. (The 
interaction results of the two links are different due to different 
arrangements of the gauges on each link. This was done to check the 
effect of gauge arrangement on the torque interaction.) Because of the 
large torque interaction effects and because the links required very close 
machining tolerances and were found to be very difficult to install and 
adjust satisfactorily, the mechanical strain intensifier was rejected as 
a practical solution to the thrust measurement problem. 
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A second attempt at mechanical intensification was made using a deflec- 
tion sensor manufactured by DSC, Inc In this device, a small cantilevered 
arm is instrumented with semi-conductor strain gauges to measure the bending 
strains produced by deflection of the end of the arm The deflection sensi- 
tivity appeared very attractive for the present application, i.e , about 
6pV/V/y-inch compared to a shaft strain of 8,85 microstrain for the 500 lb 
thrust load 

The sensor's arm has a flexure near its end to reduce the effects of 
bending moments other than those due to linear deflections It was hoped 
that this feature also would help reduce the torque and bending interactions 
in the thrust measurement Tests of several different installation arrange- 
ments showed, however, that the interactions were much worse with this sensor. 
Because of this, it was rejected as a possible thrust sensor 

3. 2 5 Direct Shaft Strain Measurements 

With the abandonment of mechanical strain intensification, a study of 
the practical problems of direct thrust strain measurement was undertaken. 

For comparison purposes, both foil and semi-conductor strain gauges were used 
In the case of the foil gauges (which have low strain sensitivity but good 
temperature characteristics), both 120 ohm and 350 ohm gauges were used 
The latter gauges were used to take advantage of the higher excitation 
voltages possible. The semi-conductor gauges (which have high strain sensi- 
tivity but poor temperature characteristics) were 1000 ohm gauges with 
gauge factors of about 155 

The gauges were mounted on the model propeller shaft and were in- 
stalled in the lathe (see Figure 3 9) and loaded in combined thrust and 
torque. The torque sensor was made up of 350-ohm foil gauges fabricated 
as special purpose torque rosettes Typical results are presented in 
Figures 3.11 and 3.12 The results are summarized in Table II For 
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comparison purposes, the expected sensitivity of the foil bridge was 
E 

.023 micro-volts/volt/lb thrust 

and the expected sensitivity of the semi-conductor bridge was 
E 

= 1 78 raicro-vo Its/volt/ lb thrust 

From the results of Table II it was concluded that the foil gauges 
are adequate for the torque measurement However, the higher sensitivity 
of the semi-conductor gauges was found necessary for reliable measurement 
of the small thrust strains Since the aircraft will produce a noisy en- 
vironment for the instrumentation system and since the sensor signals will 
be transmitted from the rotating shaft through slip rings, it is important 
that the sensor outputs be as large as possible to insure a high signal -to- 
noise ratio 

The selection of semi-conductor strain gauges for the thrust measure- 
ment was made in spite of their poor temperature sensitivity characteristics. 
The resistance of these gauges increases while the gauge factor decreases as 
the gauge temperature increases These temperature effects can be reduced 
by placing a suitable resistance (called a span resistor) in series with the 
bridge This resistor has a negligible resistance change with temperature 
As the bridge resistance increases (reducing the circuit current) , the voltage 
drop across the span resistor decreases and voltage drop across the bridge - 
increases This increase in bridge (excitation) voltage tends to compensate 
for the decrease in bridge sensitivity The value of the span resistance 
can be used to control the temperature range over which compensation is 
achieved. 

The effectiveness of temperature compensation was studied by installing 
an electric resistance heater inside the model shaft and measuring bridge 
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output as a function of bridge temperature (which was monitored by thermistors 
mounted near the gauges) . Some typical results from this study are presented 
in Figure 3 13. These results show that the uncompensated bridge has a 
temperature sensitivity of 12.3 yV/V/°F That is (for a bridge initially 
balanced at 78 °F) at a bridge temperature ofl04°F, the temperature induced 
output would be about 40% of the full-scale thrust output However 
with a 915 ohm span resistor, this bridge, at 104°F, will have a temperature 
induced output of only 5% of the full-scale thrust output. In flight, the 
propeller shaft temperature will change due to changes in both ambient 
temperature and oil temperature inside the shaft. Thus, temperature com- 
pensation of the semi-conductor bridges will be mandatory In addition, sur- 
face temperatures will be measured so that corrections in bridge output can be 
made. 

The results presented in Figure 3 14 and Table II, show that there 
IS a large torque interaction bridge output. This interaction possibly can 
be explained by one of the following arguments Theoretically there should 
be no torque strain components in the axial and lateral directions, which 
are the directions in which the thrust gauges are mounted This is true 
provided that the thrust gauges are exactly aligned on the shaft axes, and 
that the shaft is isotropic and homogenous In the present case, the ratio 
of the maximum thrust and torque strains is 1 42 4 Thus, small misalignments 
of the thrust gauges from the shaft's axial and lateral directions would cause 
a pickup of torque strain components that might be sizeable compared to the 
small thrust strains. The term misalignment here is meant to include both 
the deviation of the gauges from the axial and lateral directions caused in 
gauge mounting> as well as the deviation of the shaft axis from its theoretical 
direction. 

For the full-scale loads, the thrust and torque interactions per 



degree of gauge misalignment have been calculated to be* 


A maximum torque strain interaction of 96 8% 
of the maximum thrust strain would occur for 
each degree of misalignment of the thrust 
gauges 

A maximum thrust strain interaction of 2.23% 
of the maximum torque strain would occur for 
each degree of misalignment of the torque 
gauges 

Another possible source of this interaction could be the cross- 
sensitivity of the gauges This is normally a very small effect, but in the 
present case it may be significant because of the 1 42 4 ratio of maximum 
thrust and torque strains. 

For both the foil and semi-conductor thrust bridges, whose results 
are presented in Table II, the torque induced output was of a sign opposite 
to the thrust output. (The sense of the applied torque was the same as that 
which would occur in the aircraft } IVhen the sense of the torque was reversed, 
the sign of the torque induced strains also changed Both the foil and semi- 
conductor gauges were mounted on a shaft which had buckled slightly in an 
earlier test Before the gauges were installed, the shaft was remachined, 
but the buckling could have caused a curved shaft axis 

Another test shaft was made and a new set of semi-conductor gauges 
were mounted on it When this shaft was tested under the same conditions as the 
older shaft,' the torque interaction in the thrust bridge output was still 
present but its maximum value had been reduced to 44 3%. Also, torques 
applied in the same sense as that in the aircraft, produced torque induced 
strains of the same sign as the thrust induced strains This result is 
opposite to that observed with the older shaft This result tends to support 
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misalignment as a cause of the interaction It is very difficult to exactly 
predict the misalignment effect because four individual gauges are involved 
However, for both the old and the new shafts, the torque interaction was 
calibratable and repeatable 

A bending interaction in the thrust output was also observed. This 
interaction was much smaller than the torque interaction The bending inter- 
action can also be caused by misalignment, as well as by small differences 
in the sensitivities of the gauges making up the thrust bridges For the 
new shaft, the maximum interaction was only 7 05%. This bending interaction 
also was found to be calibratable and repeatable 

3.2 6 Instrumentation Development and In-F light Systems Description 

3 2.6.1 Introduction and Background 

A block diagram of the thrust/torque measuring system is presented 
in Figure 3 15. The elements of this system can be grouped into three sub- 
systems, Cl) strain measurement components mounted on the propeller shaft, 

(2) components located in the cockpit, and (3) auxiliary components. By 
way of background, a short discussion of the development of the strain 
measurement system is presented first 

Because of the very low thrust -induced strain levels, a high gam 
system with a high signal-to-noise ratio is required if the thrust measure- 
ment IS to be successful Initially an AC system, with a 1000 Hz carrier 
frequency, was designed and built, This system exhibited sufficient sensi- 
tivity and was quite linear However, the zero-strain bridge output could 
not be nulled, even vi?hen the bridge was resistance balanced Furthermore, 
this zero load minimum bridge output level was sufficient to saturate the first 
amplification stage. Is is thought that the problem was due to phase un- 
balance between the semi-conductor bridge arms In view of this problem, 
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the system was redesigned to use DC excitation throughout The DC system, 
which will be used for in-f light testing, has been built and is now being 
checked out on a full scale model of the propeller shaft. With the exception 
of the data recorder, all components of the instrumentation system receive 
power from a regulated ^ 15 VDC power supply operating from the 400 Hz air- 
craft electrical system 

3.2 6.2 Components on the Propeller Shaft 

i 

The thrust, torque, and bending transducers (4 -arm strain gauge bridges) 
are mounted directly on the propeller shaft The associated electronics are 
mounted on a printed circuit board that is attached to the front face of the 
slip-ring assembly (see Figure 3 16 and Section 3263 below) Bridge exci- 
tation is by a regulated, IC S VDC power supply that is part of the shaft 
instrumentation package and operates from the 15 VDC source This additional 
regulation was selected to provide a stable source for the strain transducer 
circuits 

The bridge outputs are processed by IC differential amplifiers feeding 
into operational amplifiers (See Figure 3.17). The amplifier combinations 
are adjusted to give overall gains of 1000 The amplifiers are mounted in 
high thermal resistance packages and have active internal temperature con- 
trols giving very tight temperature control This arrangement provides 
highly stable gain and low temperature drift over a range of ambient tempera- 
tures from -5S°C to + 125°C. 

3 2 6.3 Slip Ring Assembly 

A sectional view of the slip ring assembly is shown in Figure 3.18 and 
a photograph of the components, before final assembly, is presented in Figure 
3 19. The brush/ring material combination was selected on the basis of 
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information given in Reference 3 3 and after some preliminary laboratory 
checks The objective was a low resistance contact and very low slip ring 
noise The materials selected were hard brass rings and silver-graphite 
(50% silver by weight) brushes Two brushes and two internal leads are pro- 
vided for each ring to further improve the contact characteristics. 

In order to avoid the complications of a split slip ring assembly, 
the inside diameter of the assembly is 5 inches, large enough to slip over 
the flange at the forward end of the pj’opeller shaft (see Figure 3 2) A 
split spacer is clamped to the shaft, to bring the diameter up to 5 inches 
. Two aircraft quality bearings separate the inner (rotating) and outer (fixed) 
assemblies. The outer assembly is thus carried by the shaft and is not 
rigidly attached to the aircraft frame. A soft rotational restraint will 
be used, minimizing vibration effects. 

3 2.6 4 Cockpit Components 

An amplifier control box and a data recorder are located in the 
observer's (forward) cockpit The control box contains active filters, 
having 3db roll-off at about 65 Hz, and output amplifiers for use in 
adjusting the overall gam of the system (See Figure 3 20) It \tfill also 
contain auxiliary components discussed in Section 3 2.6 5 below 

The data recorder is a Gulton Industries, 8-channel, strip chart 
recorder provided by the University’s Aviation Research Laboratory (ARL) 
where it has been used for in-f light testing. It is powered by a 400 Hz 
to 60 Hz inverter, also supplied by ARL 

A bundle of shielded co-axial cables will connect the slip-ring 
brushes to the amplifier control box 

3.2 6 5 Auxiliary Instrumentation 


In order to correct the semi-conductor bridge output for temperature 
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effects, thermistors will be mounted on the shaft surface close to the strain 
gauges They will be connected to the cockpit control box through the 
slip rings and the cable Associated electronics will be contained on the 
printed circuit board on the shaft and in the control box Their output 

I 

will also be recorded on the 8-channel recorder 

The apparent thrust load will contain a contribution due to the oil 
pressure in the shaft This pressure will be monitored by a pressure trans- 
ducer in the hydraulic line between the governor and the propeller shaft. 

Its associated electronics also will be in the control box with output 
recorded on the 8-channel recorder 

Finally the control box will contain a single stage timer to provide 
a time reference for the recorder output 

3.3 Fly-by Noise Measurements 

Capability of fly-by noise measurements has been established 
Measuring equipment consists of a General Radio 1962 9601 microphone with 
associated pre-amps and power supply, feeding into a NAGRA IV-SJ tape 
recorder. Analysis of the results is accomplished by use of a General Radio 
1564-A Sound and Vibration Analyzer The noise recording and analysis 
equipment is University supplied 

Procedural checks of the noise measuring program have been satis- 
factorily completed. 

Weather conditions with sufficiently low ambient noise level to 
allow detection of the very low noise level of the Y0-3A aircraft have not 
prevailed at this writing This aspect of the program is ongoing and air- 
craft noise measurements will hopefully be achieved in the near future 
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Figure 3.1 The Y0-3A Aircraft with Cowl Removed 





Figure 3.2 The Y0-3A's Propeller Shaft and Hub 
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Figure 3.4 Approaches to the In-Flight Measurement of Propeller Performance 
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Figure 3.5 Propeller Shaft Strain Gauge Arrangements and Strain Directions 
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Figure 3.8 Mechanical Thrust Strain Intensifier, Mounted on Model Propeller Shaft 
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Figure 3.9 The Combined Loading Arrangement for the Model Shaft 
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Figure 3 13 Temperature Compensation of the Semi-conductor Thrust Bridge 
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Figure 3 14 Semi-Condtirror Thrust Rndpe Calibration for Tornno Tntnr'inti on 



Figure 3 15 Block Diagram of In-Flight Instrumentation System 
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Figure 3.16 The Instrumentation System 
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Figure 3.17 Schematic of One Typical Channel^ of On-Shaft Instrumentation 
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Figure 3.18 Sectional Sketch of Slip Ring Assembly- 
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Figure 3.19 The Slip Ring Assembly 
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PROPELLER SHAFT LOADS 
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APPENDIX A 


NUMERICAL ANALYSIS OF PROPELLERS 


by R. M. Plencner 



LIST OF SYMBOLS 


axial interference factor 
rotational interference factor 
coefficients of sine senes 
aspect ratio 
number of blades 

i 

non-dimensionai chord, C/R 
chord, feet 
lift coefficient 
drag coefficient 
diameter, feet 
induction factor 
advance ration, V/nD 
loading factor 
rotational speed, rev/sec 
torque, ft“lbs 
torque coefficient 
distance along blade, feet 
radius of propeller, feet 

vector from origin to segfhent of trailing vortex filament 
vector from origin to a blade section 
thrust, lbs 
thrust coefficient 

velocity induced at the propeller disk, ft/sec 
forward velocity of propeller, ft/sec 



normal induced velocity 

vector sum of rotational and forward velocities 
non-dimensional distance along blade, r/R 
rearward distance from the propeller, feet 

absolute geometric angle of attack 

induced angle of attack 

effective or sectional angle of attack 

blade angle 

circulation 

interval around the singularity 
efficiency 

parametric variable of the helix, radians 
tip speed ratio, V/fiR 
density, slugs/ft^ 
propeller solidity, BC/2irr 

angle between the relative velocity and the horizontal pi 
rotational speed, rad/sec 



I. INTRODUCTION 


A great deal of effort and expertise went into modeling the aero- 
dynamic characteristics of propellers in the early part of this century, 
However, with the coming of the jet age in the early - 19S0' s , the pro- 
peller was made obsolete for many applications and new theoretical work 
on propeller theories was almost non-existent Therefore, most of the 
work done in this area was undertaken before digital computers were de- 
veloped Since the problem of flow around a propeller is by nature very 
complex, these early approaches to the problem necessitated many simpli- 
fying assumptions and often made use of graphical techniques to obtain 
solutions to the problem 

Recent years have seen a renewed interest in propeller modeling for 
predicting their performance in a wide variety of applications These 
applications include more efficient propellers to reduce fuel consumption, 
propellers suitable for STOL aircraft, low noise propellers for commercial 
and general aviation aircraft, as well as rotors for wind mills. Many of 
the applications require more accurate prediction over a wider range of 
operating conditions than is possible with classical propeller theories 

The purpose of this paper is to review the classical approaches to 
the propeller problem and then dispense with as many simplifying assumptions 
as possible in order to produce a theory that more closely models the 
actual physical phenomenon 



II. CLASSICAL THEORIES 


Momentum Theory - A propeller produces thrust by increasing the 
velocity of a large quantity of air. This production of thrust is 
associated with a loss of energy which is due to the increase in kinetic 
energy^ rotational motion imparted by the torque and frictional losses 
of the propeller blades In the simplest form of the momentum theory, 
first set forth by R E Froude [1] and W. J. Rankine [2], only the 
axial component of momentum is considered; the rotational and frictional 
losses are ignored. Therefore, the theory describes the characteristics 
of in ideal propeller. 

The axial momentum theory is developed [3, 4] by replacing the 
propeller with a propeller disk that has an infinite number of blades, 
and the thrust uniformly distributed over the disk. The flow is assumed 
to be incompressible and irrotational in front of and behind the propeller 
disk. Under these conditions the change in pressure across the propeller 
disk IS equal to the total pressure head in front of the disk minus the 
total pressure head behind the disk. It can then be shown that at the 
propeller disk the increase in velocity over the free stream value is one 
half of the increase in the ultimate wake. The efficiency which is 
predicted by the momentum theory, which, is the maximum theoretical efficiency, 
is given by 


n = propeller efficiency = 


C13 


where V is the free stream velocity and v is the velocity increase at the 
propeller disk. 



Practical applications o£ the momentum theory are limited due to 
the gross simplifications that are made The theory neglects rotational 
energy losses, non uniform thrust loading losses, profile drag, blade 
interference, and compressible drag changes Its direct application, 
therefore, is only principally useful in obtaining an upper limit on 
efficiency or a crude first estimate for the performance of a propeller. 
A more general form of the momentum theory [S] includes rotational 
momentum as well as axial momentum. However, this added complexity does 
not increase its usefulness in most applications. 

Blade Element Theory - The blade element theory [5] analyzes the 
aerodynamic characteristics of the propeller by estimating the aero- 
dynamic forces experienced by each blade element as it moves through 
the air. The blade element theory is evolved by assuming that the 
aerodynamic forces acting on a given blade element are equivalent to the 
forces acting on a suitable finite wing of the same airfoil section 
moving linearly with the same relative velocity and the same geometric 
angle of attack that is experienced by the propeller blade section. The 
contributions from all the blade elements may be added up using a strip 
integration technique to get the overall propeller characteristics. From 
the geometry of Figure 1 it can be shown that the incremental thrust and 
torque of any blade element operating at a geometric angle of attack, 
a = 3 - <f), are given by 

= I BCplV^ cos <t> - Cp sin (2) 

and 

BC rpW^ sin + Cj^ cos (p) 


( 3 ) 
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2 2 2 2 

where W = V + Q r . The validity of equations (2) and (3) is largely 

dependent on the characteristics of the wing used to represent the 

propeller blade section In order to accurately account for the inter- 

ference effects, the characteristics of the assumed wing would have to 

change from station to station However, the blade element theory greatly 

simplifies the problem by assuming the wings used to represent the blade 

element sections have an aspect ratio of six. Thus the problem is 

1 

completely specified by equations (2) and (3) when the chord and blade 

angle distribution and the aerodynamic characteristics of an aspect ratio 

six wing of the same airfoil section as the blade element are known. 

The blade element theory better takes account the effect of the 

geometric shape of the propeller than does the momentum theory. However, 

it still fails to accurately account for the interference velocities 

induced by the trailing vortex system, 

\ 

Vortex Theory - The vortex theory of propellers is basically a com- 
bination of the momentum theory and the blade element theory. Therefore, 
it is often referred to as the blade element-momentum theory.' To avoid 
confusion with other theories it will be referred to as the Glauert vortex 
theory [5, 6, 7] in this paper. The Glauert theory is based on the 
assumption that the trailing vortex filaments which are produced by the 

t 

rotating blades form helical vortex sheets as they pass downstream. An 
exact application would require the induced velocity produced by the 
trailing vortex system to be computed at each blade station. However, 
Glauert simplified the problem by assuming the propeller to have an 
infinite number of blades. This assumption removes the periodicity of 
the flow, therefore, the velocity for any given radius is constant over 
the propeller disk. As a result, the momentum theory may be used to 
evaluate the interference velocities. The geometry of the problem is 
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given in Figure 2 a and a' are defined as the axial and rotational 

interference factors respectively V(1 + a) is the fon^ard velocity, V, 

plus the axial interference velocity. J^r(l - a') is the angular velocity 

of the propeller blade section minus the rotational interference velocity. 

The effective velocity, W^, of the blade element is just the vector sum 

of V(1 + a) and fir(l - a']. From the geometry given in Figure 2 the 

following relations may be defined 

> 

W sin (j) = VCl + a) (4) 

W cos $ = J^rCl - a') C5) 


It IS often convenient to define a non dimensional speed ratio, X, 
as the forward speed divided by tip speed. Using equations (4) and (5) 
the speed ratio may be written in terms of the interference factors as, 

. V 1 - a'^ 

^ ■ TTT-*“ * 

The equation for the elemental thrust and torque are found in the same 
manner as for the blade element theory, except now they are written in 
terms of the interference factors and are given in non-dimensional form as. 


R (I " (Cj^ cos (j) - Cp sin (jt) sec^ 4> 

and 

d*^c r ^ 2 2 

R ^ = a ^ (1 - a'r sin <|) + cos sec^ <p 


(7] 

( 8 ) 


where the non dimensional thrust and torque coefficients are defined by 

T = T^ ir R^ p (9) 

and 

Q = TT R^ p (10) 

respectively Equations (6), (7), and (8) are not sufficient to determine 
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the interference factors a and a'. It is now necessary to make use of 
the assumption of an infinite number of blades ^ This makes it possible to 
apply the momentum theory to obtain relations needed to evaluate a and a' . 
The results of the axial momentum theory, where rotational momentum is 
neglected, can be used to determine the axial interference factor as 


a(C,' cos (|) - Cp, sin 
a _ L D 

1 + a 4 sin Ip 

I 


( 11 ) 


Equation (IX) is not exact. In addition to neglecting the flow periodicity 
due to a finite number of blades, it also is based on the assumption of 
no slip stream contraction, and no rotational velocity component. The ro- 
tational interference factor may be evaluated by equating the torque to the 
rate of increase of angular momentum. This relation may be written in terms 
of the rotational interference factor as. 



sin 0 + cos (fi) 
4 sin <|) cos (p 


( 12 ) 


Equations (6), (7), (8), (11), and (12) are sufficient to determine 
the thrust and torque of a given propeller when the number of blades, 
the chord distribution, blade angle variation and the aerodynamic charac- 
teristics of the airfoil used at each blade location are known. 

Due to the assumptions made, the above set of equations accurately 
apply only to a lightly loaded propeller with a large number of blades. 
However, reasonable results may be obtained for many conditions varying 
greatly from the above restrictions The predicted thrust and torque 
distributions along the blade differ from the actual distributions for 
propellers operating at high blade loadings or for propellers with square 
tipped blades In the case of square tipped blades, the vortex theory 
fails to predict the fact that the thrust must fall to zero at the tip [8]. 
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Goldstein's Theory - Goldstein [9] dispenses with the assumption that 
the spacing between successive trailing vortex sheets is small and assumes 
there is a bound vortex filament springing from each blade element The 
local strength of the trailing sheet is the negative of the circulation 
gradient around the corresponding blade section The circulation must fall 
to zero at the root and the tip These trailing vortex filaments approxi- 
mately follow a helical path in the slip stream All the filaments together 
form a helicoidal surface. Goldstein specifically considers the optimum 
circulation distribution for which the energy lost due to the production of 
the trailing vortex system is a minimum, for a given thrust. Betz proved 
for lightly loaded propellers that this optimum circulation distribution 
corresponds to the requirement that the vortex sheet move rearward as if it 
were a rigid helicoid. Thus, the flow in the wake between the helicoidal 
surfaces is that of an inviscid, continuous, irrotational fluid, with zero 
circulation The circulation around a given blade element is equal to the 
discontinuity in the velocity potential of the respective point on the 
helicoidal vortex surface Since the problem is one of potential flow, it 
IS possible to apply Laplace's equation, 

(}) = 0 C13) 

where cj) is the velocity potential. Using the relations 

c = e - ^ (14) 

and 

n r 


V 


V 


(15) 



Laplace's equation in polar coordinates may be transformed to a more 
convenient form given by, 
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U 


2 


3y 


2 . . Cl - 


^ =0 

3? 


(16) 


The boundary condition for the problem becomes. 


3<j) _ _ fl^r^ wV 


( 17 ) 


for ? = 0 or 7T and 0 _<r _<R. In addition, <j) must be continuous every- 
where except at the helxcoidal vortex surface, and its derivative must 
vanish when r is infinite. Goldstein goes through a rigorous solution 
of this problem to obtain the loading factor, k, which is just a non- 
dimensional circulation defined by 


k = 


r B 

2w 7T r tan (j) 


(18) 


Lock [10] applied Goldstein's solution to a method for predicting the 
characteristics of an arbitrary propeller. The development of this method 
parallels that of the Glauert vortex theory quite closely, except that 
Goldstein's solutions are used to determine the interference velocities 
at the blades The equations that are obtained for the interference 
factors are 


1 + a 


2, O' C, cos 4 - sin tj) 
cos (j> L D 

k . 2 , 

4 sin (p 


(19) 


and 


, 2 , 0 C,, cos di + C, sin (b 

a' _ cos (b D L 


1 - a' 


4 sin (p cos (f> 


( 20 ) 


bquations (19) and (20) differ from equations (7) and (8) of the Glauert 

.2. 

PQC (h 

vortex theory by the factor — ^ It should be noted that in the limit 
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when the number of blades goes to infinity. 

k = cos^ (f) (21) 

Thus, equations (19) and (20) reduce to the Glauert vortex theory, 
equations (7) and (8), for the case of an infinite number of blades. 

Goldstein's loading factor, k, only applies exactly when there is 
a particular distribution of circulation along the blade to give optimum 
loading as defined in Goldstein's original problem. However, Lock shows 
that equations (19) and (20) make adjustments for tip losses and are 
generally more accurate than the Glauert vortex theory formulae, especially 
in the case of a square tipped blade. 

Theodorsen' 5 Theory - Theodorsen [11, 12, 13, 14, 15, 16] attacks the 
same problem of the optimum propeller as was developed by Goldstein. How- 
ever, Theodorsen expanded and extended the results of Goldstein's work 
Betz's requirement that the trailing vortex sheet move rean^ard as a rigid 
vortex sheet for an optimum loaded propeller was proved only for a lightly 
loaded propeller Theodorsen extended the interpretation of these results 
to apply to heavy loadings . 

Goldstein's original solution also applied only to infinitely light 
loadings Theodorsen showed that it is possible to extend these results 
to heavy loadings by including the rearjvard displacement velocity of the 
helicoidal vortex surface and doing all the calculations in the wake 
infinitely far behind the propeller This does not change the potential 
problem that must be solved. However, it does change the operating con- 
ditions which a given solution corresponds to. Including the rearward 
displacement velocity of the vortex surface, instead of assuming it to be 
negligibly small, takes account of the contraction of the slip stream. 



Theodorsen's method depends on the solution for the propeller loading 
coefficient, k, which is a function of the distribution of circulation 
along the blade. Goldstein's results may be used to evaluate the loading 
coefficient. However, Goldstein was only able to obtain a small nimiber 
of solutions. Iherefore, Iheodorsen used an electrical analogy method 
to experimentally measure the loading function By using the electrical 
analogy, he was able to obtain the loading function for dual rotating 
propellers as well as for single propellers. The drawback of this method 
IS that one is limited to the charts and tables which Theodorsen gives 
for the loading function. 

Application of the Prandtl Theory - Lock [17, 18, 19, 20] put forth 
a theory which solved the propeller problem through an application of the 
Prandtl theory [21] for a monoplane wing In calculating the induced 
velocity along the blade, the actual helicoidal trailing vortex sheets are 
replaced with a system of vortices in two-dimensions . This method 
assumes that the flow around the helicoidal vortex sheets can be approxi- 
mated by the flow around a system of parallel two-dimensional vortex 
sheets whose strength is equal to the actual strength of the helical 
vortices at the same radial distance from the axis of rotation. The 
distance between successive sheets is equal to the pitch of the airscrew 
The Fourier series may then be used to represent any arbitrary distribution 
of vorticity with radius Strict application of this theory would require 
the coefficients of the Fourier series to be found in a manner analogous to 
that employed in w'lng theory, where a system of simultaneous equations in 
the coefficients must be solved. However, to simplify the technique, a 
graphical method is applied to determine the coefficients of the series and 
corrections are then applied to account for errors introduced by replacing 



the helicoid vortex sheet by a system of two-dimensional plane laminae 
The method in its final form is strictly a graphical technique, using charts 
developed by Lock 

Moriya Theory - Moriya [22] developed a calculation method for the 
aerodynamic characteristics of a propeller in which he obtained the induced 
velocity by introducing an induction factor, I To obtain the induction 
factor It is assumed that the trailing vortices springing from each blade 
form a helicoidal surface which extends rearward from the blade to infinity 
The downwash velocity at each blade element due to the whole vortex system 
IS then calculated by application of the Biot-Savart law The downwash 
integral becomes singular when calculating the velocity induced at a point 
where the trailing vortex filament springs from that same point. If 
the vortex filament were straight instead of helical, it ivould exhibit the 

same singular behavior Therefore, Moriya introduced an induction factor, 

I, which defined as. 


I = 


d w 
n 

d w 


"1 


( 22 ) 


where d is the normal velocity induced by a helical vortex filament and 

d w IS the normal velocity induced from a straight vortex filament At 
^1 

the point where the vortex filament leaves the blade, the induction 
factor takes on the value of unity Induction factors for various advance 
ratios were calculated graphically by Moriya and are tabulated in reference 
[ 22 ] 

A successive approximation scheme to determine the aerodynamic char- 
acteristics of an arbitrary propeller was developed by Moriya In this 
scheme the circulation around each blade element is calculated corresponding 
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to an equivalent monoplane wing of aspect ratio six Using this calculated 
circulation distribution and the suitable induction factors, the normal 
velocity at each station is obtained by graphically integrating the following. 


w = 
n 



(23) 


An equivalent aspect ratio for each blade section can then be calculated 
from 




(24) 


This scheme is repeated by calculating the circulation distribution 
corresponding to the new equivalent aspect ratio at each blade element 
This Iteration is repeated until the equivalent aspect ratios of two 
successive approximations are sufficiently close The lift and drag coef- 
ficients are found in the final iteration and are used in the normal strip 
theory formulas, equations (2) and (3), to find the thrust and torque 
coefficients It is reported by Moriya that only one or two iterations are 
needed to obtain a good degree of accuracy 
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III THEORETICAL INVESTIGATION 


Most of the preceding theories involve many simplifying assumptions 
to obtain the interference velocities needed to analyze the aerodynamic 
characteristics of a propeller The equations and method needed to solve 
for the induced velocities with a minimum of simplifying assumptions in 
a manner analogous to that of finite wing theory will now be developed 
. The basic assumption is that of a helicoidal trailing vortex sheet which 
extends rearward from the blade to infinity. The velocity induced by a 
segment, d£, of a trailing vortex filament can be calculated at each blade 
section by means of the Biot-Savart law which is given as. 


6 V (r) = 


r d£ X Ct-s) 
4 IT 


(25) 


|t-sl^ 

where t is the vector from the origin to the blade section at which the 
induced velocity is to be calculated and s is the vector from the origin 
to the- vortex filament segment d^. For the geometry depicted in Figure 3 
the following relations can be writ’ten. 


t ~ R x"*! . (26) 

f /\ ^ 

T = R X cos 6 1 + Rx sin 9 } + Ra tan R (27) 

dZ - R X sec d0 [- cos (|) sinBi + cos <}) cos 83 + sin (j> k] (28) 

The induced velocity at a point x^due to a semi-infinite vortex filament 
springing from the point x on the blade is found by substituting equations 
(26), (27) and (28) into equation (25) and then integrating this equation 
over 0 from zero to infinity The normal component of this induced velocity 


is given as. 


dw = 
n 


r 

47rR 


/. 


Ir 2 

x|x 


xx'^ cos 91 + [x(9 sin 9 + cos 9) - x’*l y 


[X 


^2 


2 xx' cos 0 + 


x' . 


0' 


^ 2 ^ 3/2 ^ ' 2 ^ 1/2 


d 0 


(29) 
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The normal induced velocity at the point x'' due to the vortex filaments 
springing from the radius x on all the blades is 


dr 


B 


[x^ - xx'' cos 0^]x^+ [x(0 sin + cos 0^)-x''] 


"“■> ■ h I [x'2 - 2XX- cos e, ♦ x2x2)2/2 , ,2]l/2 


d0 (30) 


, Q D 2tt (B-k) 

\vhere 0, = 0 + ^ — 

k B 


The normal induced velocity at some point x" due to all the vortex filaments 
springing from each of the blades is found by integrating equation (30) over 
the blade and can be expressed as 




dr . B 

dx 


4ttR 


[X 


^-xx'* cos 0j^]x^+ [x(0 sin 0^ + cos 


- 2xx'* COS0, 


^ 2 ^ o2,2,3/2ri2^ 2,1/2 

+ X + 9 A ] [a +X-'J 


d0dx 


(31) 


Equation (31) is quite complex and is difficult to handle since the inner- 
most integral contains a singularity at x = x^ that is not shown explicitly. 
This singularity can be made explicit by applying the induction factor, I, 
introduced by Moriya [22] and defined by equation (22) The induction 
factor IS simply the ratio of the normal induced velocity for a helicoidal 
trailing vortex system to the normal induced velocity for a straight trailing 
vortex system of the same strength The induced velocity has a component 
from the vortex sheet springing from each blade The induced velocity for 
the helicoidal vortex sheet springing from the blade at which the calcula- 
tions are being made becomes infinite at x = x^ The induced velocity for 
a straight vortex system also becomes infinite at x = x" The infinite 
, velocity in both cases is a result of that portion of the vortex filament 
in the neighborhood of the point x' . Since the property of the infinity at 
x = x'' IS the same in both cases, the ratio of the two induced velocities 
must be unity at the point x = x". The contributions to the induction 
factor at X = x'' from all other helicoidal trailing vortices springing from 
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all other blades must be zero since only the blade at which the induced 
velocity IS being calculated has the property o£ infinity at x = There- 

fore the induction factor takes on the value of unity at x = x^. 

The induction factor, which is a continuous function of X, x and x' 

\ 

only, may be evaluated numerically by computer. IVhen x approaches x" the 
induction factor becomes hard to numerically determine with accuracy. 
However, since the induction factor is a continuous function and equals 
unity when x = x"*, the values when x approaches x' may be interpolated 
to a high degree of accuracy Typical curves of the induction factor 
for X = 0.3 are shown in Figure 4 

The normal induced velocity at a point x' may now be expressed 
compactly as. 


w (x"3 = 


1 dT 
dx I 


-:r dx 


C32) 


4tiR x-x 

The singularity at x = x" is now shown explicitly in equation (32) 

The method is now formulated by applying the fundamental equation of 
finite wing theory. 


a = a - a 

g O 1' 

a is the geometric angle of attack and is defined by 

a = 3 - Arctan (^) 
g X ^ 


C33) 


(34) 


The effective, or sectional angle of attack, a^, is evaluated by applying 


the Kutta-Joukowski theorem and is evaluated as. 


2 r 


R 5) (X^ + x"^)^-^^ a -C ‘C 
^ ■’ o 


The induced angle of attack, is given by 


w w 
n • n 

\ IT = W 


(35) 


(36) 
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Putting equation (343, (35], and (36) into equation (33) results in the 
fundamental integral equation which must be solved for the circulation, F 
The boundary condition for the circulation requires that, 

F(0) = F(l) = 0 (37) 

Therefore, we may express the circulation as a sine series given by 
00 

r(x') = E A sin mrx' (38) 

n=l 

Using this series, the fundamental integral equation can be expressed as 


2EA sin n nx' 
n 


Z n ttA 


n 


a c R n (X^ + 
0 ^ 


4tt R^Q. (X^ + 


✓ 

rA 

cos n 


X “ X 


3 - Arctan (— ^) 


(39) 


Difficulty now arises in determining the second term in equation (39) since 
It contains a singularity in the integrand at x = x"* This problem is 
handled by breaking the integral into three separate integrals. 


rl 6 ^x"+6 

J M dx = I M dx + I M dx + I M dx 
0 *^o x'^-S ♦^x'*-6 


(40) 


where M = 


0 

cos nTTX 


I and 6 is a small but finite value The first and 


X - X 

third integrals on the right-hand side* of equation (40) may be evaluated 

using established numerical techniques The second term, however, requires 

special consideration It can be shown that the integrand, M, will go to zero 

-xi6 


in the limit as 6 goes to zero Thus, 


f*A * ^ 

. x-d 


M dx will have a finite value. 


A method to evaluate this integral for a finite 6 is required This is 
accomplished by first making a variable substitution using, 

- X - x" (41) 


T 
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The integral can then be written as 

,6 


L 


cos Ti TT (T t X~*) 


I dx 


(42) 


Equation (42) is then expanded in a power series in x It should be noted 
that the induction factor, I, is a summation of integrals. The induction 
factor IS a continuous function by virtue of the manner in which it was de- 
fined However, each of the integrals making up the induction factor is 
not continuous and these integrals do in fact contain singularites There- 
fore, It is necessary to expand the induction factor in a power series in a 
manner which will eliminate these singularities. Since the singularities 
occur for the helix angle, 9, equal to zero, they may be eliminated by in- 
tegrating the singular integrals from 8 equal e to “> and then separately ex- 
panding the integral from 6 equal 0 to £, which contains the singularity. 

Once the power series in x is obtained, it is integrated term by term and 
2 

terms of order x and larger are dropped This gives the folloitfing approxi- 


mate relation. 


x+6 


cos n TTX 


x-d 


I dx = 26 cos (nTTx'') ^ 

k=l J 


2 rx"+ 6 

M dx + 26 


X'- 6 


x^cos (niTx^) 


^2. 

2(4 +x ) 


(1 + Zn(2c^(\^ + x"*^)) - mr 2 sin nfrx' 


5£n (6) 


4x^cos (mrx'') 


2 2 
2(x" + 


(43) 


where e is a small value and e >> 6. The singular integral can thus be 
approximated for the finite region of width 6 around the singularity. 

In order to evaluate the coefficients, A^, of the infinite sine senes 
given in Equation (38), the series is truncated at N terms. The integral 
Equation (39) is then applied at N locations along the propeller blade. 

This gives a system of N equations in the N unknowns This system of 
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equations is solved for the coefficients, A^, which are used to evaluate 
the circulation, using Equation (38) The effective angle of attack can 
then be found using Equation (35) . Once the effective angle of attack is 
found at each designated position along the blade, the lift and drag co- 
efficients can be found directly from the two-dimensional aerodynamic 
characteristics of the airfoil used at that position on the propeller 
blade The differential thrust and torque coefficients can be found by 
applying Equations (2) and (3) . 

The preceding method has been programmed for use on a digital com- 
puter. An explanation of the program is given in Appendix A and a list- 
ing of the program is given in Appendix B. 
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APPENDIX A 

Computer Program Explanation 

The computer program listed in Appendix B was developed following 
the method proposed in Section III of this report The variables in 
Table V must be inputs to the program. 

Table V Input Variables to the Computer Program 
AO - Lift curve slope per radian. 

BT - Blade pitch angle in degrees which must be specified at N locations 
along the blade. 

C - Non-dimensional chord which must be specified at N locations along 
the blade. 

CDO, CDl, CD2 - Drag coefficients of the airfoil section where the total 
drag coeffieient, CD = CDO + CDl a + CD2 a^. 

DELI - Interval around the singularity which is not integrated directly. 

EL - Tip speed ratio. 

N - Number of terms in the sine series 
NB - Number of blades. 

R - Radius of the propeller in feet. 

SKPIl, SKPI2 - The upper and lower limits, respectively, around each of 
the N blade locations, between which the induction factor 
IS interpolated 

WO - Rotational component of tip speed of the propeller in feet per second. 
X,W - Languerre-Gauss integration coefficients. 

ZTP - N locations along the blade at which the calculations. for the sine 


series are made 
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The program uses a 44 point Languerr e-Gauss integration scheme [25] to do 
the integration from zero to infinity in order to evaluate the induction 
factors. Induction factors near x’ are interpolated using a spline inter- 
polation formula [26l 



where the C's are the spline coefficients generated by the spline inter- 
polation routine developed in reference [27] . The spline integration 
technique has been found to be very accurate since it fits a curve with 
continuous derivatives through second order to the given function and then 
integrates the area under this curve. This integration scheme is also used 

a 

in calculating the total thrust coefficients from the differential thrust 
coefficients . 
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to become less accurate due to the integrand M becoming very large. A 
value of 6 between 10”^ and 10~® has been found to give good results. 

The matrix equation for the coefficients of the sine series is solved 
using Gauss elimination. The answer is checked for accuracy by substituting 
it back into the matrix equation. 



O'' s 
^ 0 

c? 

F s 

tg 


Oo 


A 


C THIS PR'JGPAN COMPUTES THE THRUST COEFFICIENT FUR AN ARBITRARY 

C PROPELLER DESIGN USING THE METHOD DEVELOPED IN CHAPTER THREE 

C OF THIS PAPER. 

C 

C ' 

C INPUT VARIAdLES TO THE PROGRAM ARE: 

C 

C AO - THE lift CURVE SLOPE PER RADIAN 

C BT - THE BLADE ANGLE AT EACH OF THE N LOCATIONS ALONG THE BLADE 

C CDO, COL, CDZ - SECTIONAL DRAG COEFFICIENTS 

C C - THE CHORD NQN-0 I MENS I ONALI ZED WITH THE RADIUS 

C DELL - THE INTERVAL AROUND THE SINGULARITY THaT IS INTEGRATED 

C SEPARATELY 

C EL - THE TIP SPEED RATIO = V/WG*R 

C EP - A SMALL NUMBER WITH EP»DELi 

C N - THE NUMBER OF POINTS IN THE SINE S-ERIES 

C N8 - THE NUMBER OF PROPELLER BLADES 

C R - THE RADIUS OF THE PROPELLER IN FEET 

C SKPIL,SKPI2 - BLADEWISE LOCATIONS BETWEEN WHICH THE INDUCTION 

C FACTOR *IS INTERPOLATED 

C WO - THb ROTATIONAL SPEED IN RADIANS/SEC 

C X,W - LANUUERRE-GAUSS INTEGRATION POINTS 

C ZTP - THE N NON-DIMENSIONAL SLADEWISE LOCATIOivS AT WHICH 

C THc CALCULATIONS ARE MADE FOR THE SINE SERIES 

C 

IMPLICIT REALMS (A-H,0-Z) 

INTEGER-4 IDUMI9 J 

DIMENSION ZTP(8L), ZT(8l), XZ(3I), Qt81), AA { 4 , 8 L ) , SS I 5 , 3 L) , 

1 QINT(9 ),SINT(9 ),S(ol), RINTI9 ),ALP019 },C{9 ), 

I YH81), CA(9 ,9 ), RHS(9 ), ASI9 ), CCA(9 ,9 ), 

3 15(315, TIGRL(U,LL), DTL(1L5, BETA(115 

CQMM0N/BBB/X(44) , W(445 
CQMMON/CCC/ TT, TTP , PI, EL, EP, ND 
C RADIANS TO DEGREES CONVERSION 

R0,C = 57.2957795L3082300 
PI =3.1415926535897900 

C READ IN DATA FOR A 44 POINT L ANGUERRE-GAUS S INTEGRATION 

READ(5,5) (X(L), W(L), L=I,44) 

5 FORMAT (2D30.0) 
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WRITE (6,6) (X(L), W( L) , L=L,44) 

6 FORMAT (• », 030-15, 030.15) 

N=9 

NB = 3 

EL=0. 293600 
EP=0-30D0 
W0=153. 93300 
A0=O.07D0-ROC 
R=1.40D0 

000=0.0333432100 
C01=~0 .O0SOO893O0*R0C 
C0 2=0.000633’93DO>!=RDC*=1'2 

C NP IS THE NUMBER OF POINTS USED IN THE SPLINE INTEGRATION 

C NP MUST BE AT LEAST 4*N SO THAT WHEN INTEGRATING FRON X=0 TO 

C ZTP-DEL OR ZTP+DEL TO 1 THERE WILL ALWAYS BE AT LEAST 4 POINTS 

C FOR THE SPLINE INTEGRATION. (SUBROUTINE SLPIZ WILL FAIL WITH 

C LESS THAN FOUR POINTS/ 

NP=(N-1)=«^10+L 
00 9 J = 1 ,N 

C READ IN. THE BLADE ANGLE, CHORD, AND CORRESPONDING LOCATIONS 

C ON THE BLADE 

RE AO (5, 7) BT, C(J), ZTP(J) 

7 FORMAT (3F10.6) 

BETA! J )=BT/RDC 

C ALPG IS THE GEOMETRIC ANGLE OF ATTACK 

IF ( ZTP ( J ). L E.O .000100 ) AL PO ( J ) =' B T/ RD C-P I /2 .0 00 

IF I ZTP { J ) . GT .0 .000 IDO ) ALPO(J)= B T/ROC-DATAN( EL/ZTP ( J ) ) 

WRITE 16,3) ZTP(J), BT, ALPG(J), C(J) 

8 FORMAT ( » • , 4F20.6) 

9 CONTINUE 

DO LO 1=1, NP 

10 ZT ( I } = FLOATi I-n/FLOAT(NP-i) 

0EL1=0 .OOOOLIDO 

DO 400 J=1,N 
DO 4 I=1,NP 
IS I I ) = 0 
4 YI I I ) = 0, 000 

WRITE (6,11) ZTP(J) 

11 FORMAT (*0*, »2TP=», F8.3J 
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SH=EL**2+ZTP{ J)**2 


THIS SECTION CALCULATES ALL THE INDUCTION FACTORS NEEDED 
TTP = ZTP( J) 

tt=ttp 
N0PT=2 

PL2 IS THE SUM OF THE INDUCTION FACTORS FOR 6LA0ES I AND 2 
DIVIDED BY (ZT-ZTP) 

CALL CALCI (P12, NOPT) 

NOPT=3 

P3 IS THE INTEGRAL FROM EP TO, INFINITY OF THE INTEGRAND M 
rt'HERE I3=( ZT-ZTP ^INTEGRAL OF M 
CALL CALCUP3, NOPT) 

NOPT-L 
TTP=ZTP(J) 

NAB = 0 
NCTR=0 
ICTR=0 

C INTERPOL.ATE VALUES OF I BETWEEN SKPIl C SKPI2 

C WHERE SKPU fl SKPI2 ARE READ IN FOR EACH ZTP LOCATION 

READ[5,12) SKPIl, SKPI2 

12 FORMAT (2F10.3) 

WRITE SKPIl, SKPI2 

13 FORMAT {'0*,' SKIP DIRECT I CALCULATION BETWEEN ZT=', F7.3, 

I • AND' , F7.3 ) 

DO 15 I=1,NP 
TT = ZTm 

C find THE VALUE OF THE DO LOOP PARAMETER FOR WHICH ZT=ZTP 

IF {DABS(TTP-TT) .LT.O.OODOOIDO) ICTR=I 
IF (TT.LT-SKPI1.0R.TT.GT.SK°I2) CALL CALCKXI, NOPT) 

C STORE THE INDUCTION FACTOR IN THE ARRAY YI 

IF ( TT.LT.SKPIl.OR.TT.GT.SKP 12 ) YI(I)=XI 
C STORE WHICH INDUCTION FACTOR CALCULATIONS HAVE BEEN SKIPPED 

IF (TT.GE.SKPI L.AND.TT,LE.SKP12 ) ISm = l 
15 CONTINUE 

DO 17 1=1, NP 
ICK=0 
TT=ZT( I ) 
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IF ( ISI I ). NE. 1. AND* TT.GE.sk PI 1-0. 0 74D0. AND. TT, LE . SKP I 2+0 . 074D0 
1 .OR.I .EQ.ICTR) ICK=2 

NAB IS THE NUMBER OF POINTS AT WHICH THE DIRECT CALCULATION 
OF THE INDUCTION FACTOR HAS SEEN CALCULATED 
IF (ICK.EQ.2) NAB=i\IA8 + l 
IF {ICK.EQ.2) XZINAB3=TT 

STORE THE COMPUTED INDUCTION FACTORS IN ARRAY Q 
IF nCK.EQ.2. AND. I.NE.ICTR) Q(NA6)^YI(n 
IF (I.Ey.ICTR) 0[NAB)=1.0DO 

17 IF (I. EQ.ICTR) NCTR=NA8 
WRITE (t>,21) NAB, ICTR, NCTR 

21 FORMAT <'0', *NAB=», 3110) 

WRITE (6,19) (XZIL), QIL), L=1,NAB) 

19 FORMAT ( ' • , 2F20 .5 ) 

USE THE SPLINE INTERPOLATION FORMULA TO OBTAIN THE INDUCTION 
FACTORS THAT WERE SKIPPED IN THE DIRECT CALCULATIONS 
SLPIZ IS A COMPUTER SUPPLIED SUBROUTINE WHICH WHICH GENERATES 
SPLINE INTERPOLATION COEFFICIENTS. 

CALL SPLIZ (XZt Q, NAB, AA, SS ) 

DO 18 I=>1,NP 
TT = ZT{ I ) 

C USE SPLINE INTERPOLATION FORMULA 

IF (IS ( I I .EQ. 1 . AND. I .LT. ICTR) 

1 YI(I)=AA( 1,NCTR-1)*(XZ(NCTR)-TT)<‘=!‘3 

1 +AA{2 ,NCTR-l)« (TT-XZ(NCTR-i) 

2 +AA(3,NCTR-U=S=(XZ{NCTR)-TT) 

3 +AA(4,NCTK-l)>i=(TT-XZ(NCTR-l) ) 

IF { I .EQ. ICTR) YI(I)=1.0D0 

18 IF t IS (I 1 . EQ. I . AND. I .GT. ICTR ) 

L Y I ( I ) = AA ( I , NCTR )*(XZ(NCTR + 1)-TT)<'’!'3 

2 +AA(2,NCTR)*^(TT-XZ(NCTR) 

3 +AA( 3,NCTR )’:'{ XZ(NCTR+L)-TT) 

'4 +AA(4,NCTR)>^ITT-XZ{NCTR) ) 

IF ( ZTPI J) .GT.DfcLL) TT=ZT P (J ) -0 L LI 
IF (ZTP( J) .GT.DELL) 

1 XI BZ = AA( 1 ,NCTR-1 )^(XZ(NCTRj-TT)4:*3 

1 +AA( 2,NCTR-L)=<= (TT-XZ( NCTR-1 ) 

2 +AA(3,NCTR-1)«(XZ(NCTR)-TT) 

3 +AA(4, NCTR-U* (TT-XZ (NCTR-I) ) 
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IF (l.OOO-ZTP(J)-GT.OeLU TT=ZTP ( J ) +DEL 1 
IF t 1.000-ZTP{ J) .GT.DELl) 

1 XIAZ=AAU ,fJCTR )’i'(XZ(NCTR+l)-TT)#-3 

2 + AA ( 2, NCTR)>:'(TT-XZCNCTR 3 

3 +AA(3,NCTR)=i‘(XZ{N!CTR+-I)-TT) 

4 +AA{4,NCTR)’f'(TT-XZ(NCTR) ) 

C ALL THE INDUCTICiN FACTORS HAVE iNO.w b£EN CALCULATED AND 

C STORED IN YI, XI6Z, XIAZ. XI3Z C X I AZ ARE THE INDUCTION 

C FACTORS AT ZTP-DEL AND ZTP+DEL RESPECTIVELY. 

DO 210 NS=L,N 
APG=NS*PI*ZTP( J ) 

WRITE (6,20) NS 

20 FORNAT (’O', *+++ NS = ', 13, ’ + + -(- + + ++ + +++•) 

DO 22 IJ=1,N 
QINT( I J)=O.ODO 
RINT( IJ)^O.ODO 
22 SINT( I J)=O.ODO 
DEL=0£L1 
ND = 0 

IF (ZTPt-J).LE.DEL ) GO TO 53 
C 

C THIS SECTION CALULATES THE INTEGRAL FROM 0 TO ZTP-DEL 

C 

DO 25 1=1, NP 
ND=ND+ 1 

IF (ZT(I )-ZTP( j) .GE.-OELI ) GO TO 30 
XZ ( I )=ZT {[ ) 

C STORE THE INTEGRAND IN ARRAY 0 

25 on ) = OCQS(NS*PI’^ZT( I ) )^YI ( I) / ( Z T (I ) -ZTP ( J ) ) 

30 XZ(NO)=ZTP( J)-DEL 

U(ND> = DCOS(NS-!<P I*XZ (ND ) ) *X 13 Z/ ( XZ( ND ) -ZTP ( J ) ) 

ND1=ND-1 

IF (NS.EQ.I) WRITE (6,35) (XZ(L), Q(L), YI(L), L=L,ND1) 

IF (NS.E(3.1) WRITE (6,35) XZ(NU), Q ( ND ) , XI3Z 

35 FORMAT ( ' ’ , 3D15.5 ) 

C INTEGRATE Q USING SPLINE INTEGRATION FORMULA 

CALL SPLIZIXZ, g, ND, AA, SS ) 

SUMG=0.0D0 
ND l = ND-i 
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DO 40 K=l,N0l 

40 SUMQ=SUMQ + ,25D0=:<AAa,'K}«tXZ{K+L }-XZ(K) ) **4 

1 +.2 5D0’!'AAl 2 ,K)*(XZ(K+1 )-XZ{K) )*^4 

2 +.50D0«AA(3,K)«(XZ{,X*-1)-XZ(K) )**2 

3 +.5000*AA(4,K)*(XZ(K+L)-XZ{K) )#=!'2 

QINT(MS)-SUM0 

WRITE (6,50) SUMQ 

50 FORMAT (• 'INTEGRAL OF Q= ' , Ell. 3) 

53 CONTINUE 

IF ( IGT.EO.l ) GO TO 75 

IF (ZTP(J).LT.DEL.OR.UODO-ZTPUKLT.OEL) DEL=DEL/ 2 .00 0 
THIS SECTION CALULATES THE INTEGRAL FROM ZTP-OEL TO ZTP+OEL 
ATERM=OCOS{ ARG ) =i'2.0D0=i«DEL*PL2 

BTERMl = 2.OD0*0EL*( IZTP( J )*DCGS{ ARG) ) / I DSu/RT ( 2. ODD ) =i=SH ) 

1 ><=(l.ODO + OLOG{2.0DO»!'EP>>=>!'2^SH) ) 

2 -NS^PI=f=DSQRH2.0D0I*DSIN{ARG) ) 

6TERM2= -0EL« OLOG(OELL)^(4-000=('ZTP{J)=!'DCOS( ARG) 

1 * /{DSQRT (2.000)«SH> ) 

CTERH = 2.0DO>^DCOS( ARG)=^DEL«P3 
RI NT(NS)=AT£RM + 8T ERM I +-BT ER M2 +CT EkM 
IF (NS.EO.l) vVRITE (6,55) P12, P3 
55 FORMAT ( ‘O' , 2E20.5) 

IF {NS.EQ.D WRITE (6,60) ATERM, 8TERM1, 3TERM2, CTERM 
60 FORM.AT ( 'O' , 4E20-5) 

WRITE (6,70) RINT(NS) 

70 FORMAT (' *, 'INTEGRAL OF R=', £11.3) 

DEL^DELl 

75 RIiNT(NS)=O.ODO 

IF ( 1.0D0-ZTP( J) . LE.DEL) GO TO 160 

THIS SECTION CALULATES THE INTEGRAL FROM ZTP+DEL TO 1,0 
ND=1 

XZ(NO)=ZTP (J ) +DEL 

S(ND)=DC05{NS’i-PI*XZIN0) )*XIAZ/{XZCND)-ZTP( J) ) 

DO 100 1=1, NP 

IF { ZT (I )-ZTP ( J ) .LE.DELl ) GO TO 100 
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NO=ND-H 
XZ(ND) =ZT( I ) 

C STORE THE INEGRAND IN ARRAY S 

s(NO)=ocus(NS#Pi + zT( n )>f=Yi ( I j/(ZT( n-zTP(jn 
100 CONTINUE 
110 CONTINUE 

IF (NS.EQ.l) WRITE (c,35) XZ { I ) r S(l), XIAZ 

LNO=NP“ND 

DO 115 L=2,ND 

LL=L+LND 

IF (NS.EQ.l) WRITE (6,33) XZ ( L) , S(U, YI(LL) 

115 CONTINUE 

C INTGRATE S USING SPLINE INTEGRATION FORMULA 

CALL SPLIZIXZ, St NO, AA , SSI 
SUNS^O.ODO 
NDL=ND-1 
DO lAO K=1,ND1 

140 SUMS=SUNS+.2 500*AA{ 1 , K ) ( X Z ( K+ 1 ) -X Z ( K) )«#4 

1 +.2500>i‘AA(2,K)*(XZ(K+i.)-XZ(K) ) *’!^4 

2 * +.50D0-XAA( 3 ,K)*(XZ( K+1 )-XZ(K) )=x«2 

3 +.50C0*AA(4,K)*(XZ(i<+l)-XZ(K) )«*2 
SINT(NS)=SUMS 

WRITE (6, 150) SUMS 

150 FORMAT (• », 'INTEGRAL OF S= ' , E12.4) 

C TINT IS THE TOTAL INTEGRAL FROM ZERO TO 1 

160 TrNT=QlNT(NS> ^-R I N T ( NS ) +S I N T I NS ) 

TIGRLI J,NS)=TINT 
WRITEt c 170) TINT 

170 FORMAT ('O', 'TOTAL INTEGRAL FROM 0 TO INFINITY=», E20.5) 
C SET UP THE MATRIX EQUATION 

CA(J,NS)=2.0D0*DSIN(NSMPI=!=ZTPIJ) ) / AO-NS’i-CI J ) ’i'T I NT/A.ODO-^ 
CCA(J,NS)=CA{ J,N5) 

210 CONTINUE __ 

RriS( J ) =W04R^^«OSQRT( SH)*ALPO(J)*ClJ) 

WRITE (6,350) 

350 FORMAT (*0* , ' 

400 CONTINUE 

WRITE(6,500) ((CA(MtL), L=1,N), M=1,N) 
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500 FORMAT (« •, 9Fil.5) 

WRITE i6,550> /RHSlLJt L=1,,'V) 

550 FORMAT (*0*, 9F11.3) 

- SOLVE THE Matrix equation for THE AS VAt iicTMr A 

CALL*^GAUSZ^^rA^^° GAUSS ELIMINATION TECHNI QUE. 

N» N, RHS, ASr IDUM, I ER ) 

WRITE {5,oOOI USIL), L=1,N) 

600 FORMAT (* 'A=', D20.5) 

WRITE (6,o25) 

625 FOR.tAT CJ', |lx, J6X, 'GAMA', 16X, "ALPHA !', 

ST IGR=0. ODO 
5GAM=0.000 

VC = hO*R^^H-OSQRT(£L**2 + ZTPlNK)«*2) 

VR IS THE RELATIVE VELOCITr 
VtA,= VC/ R 
DO 650 K=1,N 

ST I GR= STIGR + TIGRL(MK,K)*AS(KJ=)=K 
CALULATE THE CIRCULATION 

o5 0 S&AM=SGAM+AS{K)*DSIN(K^PI=*7TP(MK) I 
ALI=STIGR/(4.0D0^VR) 

CALCULATE THE EFFECTIVE ANGLE OF ATTACK 

AL0=2.000;''SGAi'V{A0’4'C(NK)vvC) 

PHI=B6TA(HK)-ALQ 

CD = CD0 4-C0I-!'ALa*CD2=^^ALQ^«=2 
CY = CL*DCOS ( PH [ )-CD*DS INC PH I) 

CALCULATE THE DIFFERENTIAL THRUST 

WRITE (6,6 75 } ^ZTP J 

675 FORMAT [• V, sLo^S) ’ DTC(NK), PHI 

700 CONTINUE 

FORMULA^^ thrust over the dLAOE USING SPLINE INTEGRATION 
CALL SPLUIZrP. DTC, M, aa , SS ; 
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NMl=iN- L 
SUMT=O.ODO 
DO 70Z K=l,NMl 

7 02 SUMT=SUMT+. 2 5000* A A {2r!0*{ZTP(K + L)-ZTP(K))**4 

1 ■»-.500D0*AA(4,K) *( ZTP(K + i)-ZTP (i<)) **2 

2 +»25000*AA{l,K)*(ZTP(K + U-2TPiK) )**4 

3 50000*AA( 3, K)*( ZTP<K^-n -ZTP <KJ ; **2 
WRITE (6,705) SUMT 

705 FORMAT ('O', ‘THRUST COE FF IC I E NT= ‘ , F20.7) 

CHECK TU SEE IF THE MATRIX EQUATION WAS SOLVED CORRECTLY 
00 725 1=1, N 
CKRHS = 0.0 
DO 710 J=1,N 

710 CKRHS=CCA( I , J ) *AS ( J) +CKRHS 
725 WRITE (6,730) CKRHS 
730 FORMAT (‘O', F20.3I 
RETURN 
END 




SUBROUTINE CALC I (XI, NOPT) 

THIS SUBROUTINE CALCULATES THE INDUCTION FACTORS AS FOLLOWS: 
NOPT = 1 CALCULATES 11+12+13=1 

NOPT =2 CALULATES THE INTEGRAL OF P1+P2 FROM 0 TO INFINITY 

CALULATES THE INTEGRAL OF P3 FROM EPSILON TO INFINITY 
WHERE I=(ZT-ZTP}*(P1+P2+P3 ) 

IMPLICIT REAL*8 (A-H,0-Z) 

COMMON/AAa/ KK, NNGPT 

COMMON/CCC/ TT, TIP, PI, £L, EP, NB 

DIMENSION AN (3) 

EXTERNAL AUX 
NN0PT=N0PT 


00 10 KK=1,3 

IF (NOPT.BQ. 2- AND.KK.E0,3) GO TO 20 
IF (NOPT . EU . 3 . AND. KK-NE . 3 J GO TO 10 
CALL GLOU (AUX, ANS ) 

10 AN(KK)=ANS 

IF (NOPT .EO.l. AND. DABS (TT-TTP) ,LE. 0.00 100) 
i >:i = l.ODO 


Z£-V 



ooooo noo 


1 


IF (NQPT.EQ. l,ANO.OABy(TT-TTP) .GT. 0.00 IDO) 

1 XI=lAN(d^)+AN(2 )-t-AN(3 ) )*(TT-TTP) 

20 IF (M0PT.EQ.2) X I=AN { I ) +AN { 2 ) 

IF (iNOPT.EQ.3) XI=AN(3) 


RETURN 

END 


SUBROUTINE GLQU IAUX,ANS) 

THIS SUBROUTINE SETS UP THE L AN CUE RRE-GAUS S INTEGRATION 
IMPLICIT REALMS (A-H,0-Z) 

CQMM0N/B8B/X{44) , W ( 44 ) 

200 ANS=O.ODO 
S=O.ODO 

DO 201 1=1,44 
Y=X( I ) 

CALL AUX(Y,Z) 

201 S=S+Z*W( I ) 

ANS=S 
RETURN ' 

END 


SUBROUTINE AUX (PHI, FX) 

THIS SUBROUTINE CALCULATES THE INTEGRAND OF PL, P2, & P3 
WHERE THE INDUCTION FACTOR = ( ZT-ZTP ) >!= I NTGRA L OF (PI, P2, CP3) 
NOTE THAT THE INTEGRAND OF P IS GIVEN BY A/B BELOW 
IMPLICIT REAL-i'8 (A-H,0-Z) 

COMI'ON/CCC/ TT, TTP , PI, EL, EP, NR 
COMMCN/AAA/ KK, NOPT I 
> A IF (N0PT.EQ.3) PHI=PHI+EP 

[rTiftr P) PH IK=PHI4-2. 0DQ*PI*{ NB-KK > / NR 

A=^T*>f=2)* ( TT-TTP*DCQS ( PHI K) ) + (TT*(PHI^DSIN( PH U ) +DC0S ( P H I K ) )-TTP) 
1 vEL**2 

6={ DSORT { TT4«2 + TTP>k^ 2-2.0DO*TT*TTP*DC0S ( PHIK)-nPHI-l‘’!‘2=^EL-«2))*«3 
U’ 1 #DSQRT( £L*>^2+t:^«2 ) 

C = A / B LT-f-fP 

FX=DEXP{ PHI )4C 
RETURN 
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APPENDIX B 


FORMER APPROACH 

TO THE PROPELLER NOISE PROBLEM 



by C. 

. J Woan 


LIST OF 

MAIN SYMBOLS 

a 

0 

speed of sound 

t 

observation time 

I 

acoustic intensity 

V 

the velocity of source 

£[r) 

thrust per unit length 


forward velocity of propeller 

M = 

V/a„ 

X 

field point 

= 


X = 
0 

: Cxqj + Xq 2 + Xqj) the distance 




between the field point and the center 

Mt = 

r fl/a 
t o 


of the propeller at t 

\ = 

M-R/R 

y 

source point 



0 

the instantaneous angle the source 

P' 

acoustic pressure 


makes with J 2 " axis 


acoustic pressure due to 

0 

n 

initial angular displacement 


thrust 

u 

N 



'j' 

tan'F = ^02^^01 

R = 

X - Y 

T 

retarded time 


radius of propeller 

Q. 

angular velocity of propeller 

T 

total thrust 





[f(t)] E f(T) 



Source Point 



Field Point 


Figure 1 Coordinate System 
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2. Acoustic Pressure Perturbation 


Following Ffowcs Williams and Hawkings, the acoustic pressure perturbation 
at the field point is given by 

P. 


P =- 


I 3 


4Tl j f R.| I - MkI 


] dS 


th 


where dS is the blade surface area and P^ is the i component of the 
force acting on air by the blade surface 


3 Noise Due To Thrust 

Let -t(r) be the spanwise distribution of the thrust along the lifting line 
per unit length Then 




p'=- 


I 3 




at ' i-Mr. R.‘3tjJ 




r — ^ 


R c I' Mpy 


Pjy-.-'iO l-H 

R I- 




dr 


where P^ = -t(r) 


P2= 0 


0 


Therefore the pressure perturbation due to thrust is 


Pt 


where 


J ft ( s.) ^ C d r 

t 0 


it ’ 


(X,o 

TV =-2T 


P- 


/-> 4- t- ^ Xo 




9Cc 


A. 
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lev*) =-^ =l*(r*) r*n- rt) 

T IS the total thrust 
and 


ft ^ = Mt ^ Cf'l" 3 I Sint cost— 

[KJtl-Mj" 


"it 

+ r M-t >1!, Stri'^^ 


cosce] 




3 + r*Mt Mj si'i't — 


- ')^<, «s f — ^ 




1 


where the angle 0 is uniformly distributed over (0,2vl Then the proba- 
bility density function for the random variable 6^ is 


■f(Oo) = 


I 


O e < ZX 


0 


otherwise 


Therefore, the mean square of P 


t 

t 


IS 


Ecn > 



The intensity 
T ==: 


Then define 


liKl 

Po 


4(- I 

1 - -j-2 


Fcr, r*) fer"") Jer, ) ar dr/ 


£> 0 



where 


27t 

= -^3 j ft(rt9.) ftCn'J e.) d0„ 


The total thrust is given by 


T = if i?cr ) r c I- r ) dr 


Letting 


T - 1 = [ i\r^) dr 


CD 


the necessary condition that I* be an extremum, subject to the constraint 
J = const (i.e. constant total thrust) is 


TCrr, ) i 07 )dr; =--5- r O- r ) 

-^0 

\ifhere b is a Lagrangian multiplier 

and Z* Cr*) and b can be obtained by solving equations (1) and (2) 
Numerical Method 

The equations needed to solve this problem are 


I - ^ 


r") dr 


f- 


(B) 


(C) 


CD)' 


j'FCrr/;')J(r)dn* =-#r^O-r; 

“^0 

^ 27L 

r,*) = TI^J ^t.(rU.}-^(r,%jde, 




CislO] 


C2) 




, cosCe) ^ ^ sincej -x-o cos-^/^ 
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2 . 2. 

(E). t -t '^t Sin'f + r^Mt 

-X- 

-zy-t r^Mt cos c Op + 

CF) re] 0, t^- 

CG) * [ 1 - MpJ = ( I - h| ) - cos -I' ~ r^Mt s^n-Y SmCOj) 


The solution proceeds as follo^^s 
Given Mp, x*, and T, find b and 


1 The integrals in equations (A) , (B) , and (C) are calculated 
by Gaussian integration (Gauss -Legendre formulas or Gauss- 
Chebyshev formulas] . 


Equation (E) is solved by the Newton Method 

-V M it 




n. 


crj„- 


■fCirt,.-,) 




^ y. 


- 2Tot e.n'^ cosCOo -ft CRJ) J 

After obtaining [R^^] from (2], the value of [R*] is substituted 

into equations (A) , (B) , (C] , (D) , (F] , and (G) 


4. After applying Gaussian integration to equation (B) and assign- 
ing n different values for r* (o^rgl, n = the number of points 
used in Gaussian integration), equation (B) can be reduced to 

a system of n equations with unknowns Z* (rp (i=l,2,--,n) 

til 

[r* are the zeros of the n^ polynomial used in the Gaussian 
integration] 
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In matrix form 
AL = C 

A, nxn matrix 
C * column matrix 

5 The Lagrangian multiplier is obtained by substituting £* (r 
obtained £tom (4) into equation . 


^ H 
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APPENDIX C 

PRESENT APPROACH TO THE PROPELLER NOISE PROBLEM 
by C. J. Woan 


LIST OF 

A Fourier coefficients of F 
1 

speed of sound 

B number of propeller blades 

C complex magnitude of the 
th , 

m harmonic 

dD drag per unit length 

dL lift per unit length 

I acoustic intensity 

J Bessel function of order m 
m(Z) 

and argument Z 
M = V/a 

0 

= M • R/R 

P-P^ acoustic pressure 

^rms ^o°t-mean- square pressure 

R = X - Y 

n 2 2 2 

R = X +y +z 
o ^ 

r^ radius of propeller 


SYMBOLS 

T total thrust 

t 

t observation time 

V the velocity of source 

Vp forward velocity of propeller 

induced velocity 
X field point 

Y source point 

induced angle of attack 

(j) initial angular displacement 

0 the instantaneous angle the 

source makes with the y^ - axis 

fluid density 

r spanwise circulation distribution 

X retarded time 

angular velocity of propeller 
[f(t)J = .f(T) 
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2. Sound Pressure Due to Moving Lifting Line 

The force acting on the air at Y = J^x+r, by the moving lifting line, is 

developed as follows 

Consider a blade element at location r 



(dL COS Y - dD 71 't) , dpg == clLsiTvYd-dDcoST^ 


The force, dF, acting on the air at Y = flx + r is 
df = 6f^ Z -dpg S)TV6 J + dFg cos & 


The total sound pressure generated by the moving lifting line is given 
as 


CP-Po)(x,i) 


r r 

L Mp.) 


] 


3 Periodicity of the Acoustic Pressure Disturbance 

Consider a propeller whose hub is moving with constant velocity Vp 
Choose fixed axes 0 such that the origin coincides with the hub 

(i) at time x = 0 

(ii) the X axis is parallel to the axis of rotation 
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(ill) the hub velocity vector lies in the plane and in x-direction. 

Let the field point x~ be at (x,yjZ) and a typical source point be situated 
at Y(t) This can be expressed as the vector sum of the hub location and 
the rotational position 


Y<r) = Vf-J: + r iV , nz) = ( o, rcoss, rs,a0). 

Here 6 is the instantaneous angle the source makes with y axis This in- 
creases linearly with time, and can be expressed as fix+<j) where Q is the 
angular velocity of the source, and 41 is the initial angular displacement. 


z ^ -t - 


CRJ 

CLo 


If attention is confined to source times near t= 0 , [R] can be approximated 

CRJ (x-Vp)^ + ( ^ - rcos[0j)^ i- ( ? - rsia[G3)^| 


R-o C I + ^ ) 




r 1 1 1 ^ 

where P.o ~ 


and 


a = 1^ - 2rC'^coscej + :S SmC&j + x 


Therefore, 


[P-J = C + 2s,aC0J) 


jnt - sx-p =ilrc / -Mor) - Mr C^cojM + 


and 
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where 


H 


or 


X Vp 
Q-o R.0 


and 


H, 


sir 

(Xn 


This equation shows that as 9 increases by 2 it, fit only increases by 2 it 
Hence the basic observed frequency is w=fi/ 

and therefore has a period 27r/o) 


4 Fourier Analysis 


2ir 


(P-P^) (^,t) has a period — , where « = fi/ (r-M^^) 


and M 


V_x 

F 


0^ ^0^0 


th 


Then, defining the complex matnitude of the m harmonic in the usual way 


Cwi ( X ) — + <-" b -m ~ 


CO 

7t 


CO 


r ^ _ cm loi ,, 

\ (P-P.)cx,t)e dt 

•Jf) 


=-I_ C ^ T 

4-iL ) 


vjhe re 


I- = ~ 

7C 


n 






. iiacot , 

J e dt 


changing variables back to the retarded time, x, gives 


m 


~T” . - ffO 

~ 7C 


j j ) dr 


u;Kere dt = [ l~ dZ 

= R .0 - cos C0J +2 s.acaT) 

Since the far field approximation is in force, the only significant 
dependence of upon is through the term in the exponent. « 
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Therefore 


rrpy 

■ _ -t >n r ^ p a 1^ 0 4 ntco(T+-^) dC 

a„7t: J 99c^ J ^ ° /!^ 


-t tn CO £ 


{ m to < c^ 


C -v I* » 7 

m “ "~7 ~ - 

2 JC «. /2«(I-M„) 


( 1 -’- 


I- Kr JL 
M-h 


-r / Mi p . 

^^0 . 


Extension to B Blades 

For an assembly of B identical blades, the usual phase arguments show that 
only harmonics which are multiples of B survive, and all the rest cancel 
Hence only harmonics of the propeller blade passing frequency, Bco are 
present in the acoustic field of the complete propeller, and the magnitude 


of the m harmonic is EC 

mg 


Cm - 


I 'm & CO & 

Zit ao do (l-Mor) 








Then, the r m.s pressure is given as 


T V?-. 


The total intensity at 1C is, 


1 
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and the total thrust is 


Tt =5 (r^ + P^p dr 

0 

Method of Solution 

The problem to solved can be stated as follows 

Find the circulation distribution along the lifting line such that I 

minimum and constant, 
t 

The solution then proceeds as follows. 

Expand F as Fourier series with coefficients then 

i ~ I C , /As • ’ 'J) = I At } 

Tt " Tt C Ai 5 A 2 , A 3 • * 0 “ Tt ( A t) 


Let J = I - b(T^ - \ 

where b is a Lagrangian multiplier. 

A necessary condition that I be extremum is 


= 0 


= 1 , 2, 3, 


Solve equation CA) and (B) for and b 



